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ArP  4Hl7<-f  w 

Srruta  f«r  tli<«  Report  WlS-ONR-12  **The  Optimum  Free 
Theory  of  I/iqulds"  • 15  October  1954. 

(1)  In  e^ueiion  {7)  end  (3)  replace  all  py  by 

(2)  Th.  «».ctu..C^V(«iy.yj)  -*V(s;y,yj)l  occurring 

in  eqoetlon.  (13)t  (15),  (16),  (21),  and  (22)  should  be 
replaced  by[  ^r(Wyj)-25*rC3',  V.Vj)]. 

(3)  The  factor  multiplying  the  bracket  in  equation  (13)  should 

be 

(4)  The  factor  multiplying  the  integrals  in  equations  (15)  and 

(21)  should  be  tT  7 

(5)  The  factor  xiKtltiplying  the  integrals  in  equation  '{l6)  and 

(22)  should  be  Zn^. 
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THE  OPTIMUM  FREE  VOLUME  THEORY  OF  LIQUIDS 

by 

§ 

John  S.  Oahler  and  Joseph  O.  Hirschfelder 
ABSTRACT 

An  iterative  method  for  solving  Kirkwood's  integral  equation  fur 
the  free-volume  of  a liquid,  is  described.  This  procedure  uses  the 
results  of  the  well-known  Leonard -Jones  and  Devonshire  theory  of 
liquids  as  its  starting  point.  The  results  of  the  calculations  sug- 
gested would  provide  a good  test  of  the  validity  of  the  cell  theory 
of  liquids.  It  is  belie  vad  that  this  "optimvun  free -volume  theory" 
will  remove  several  of  the  undesirable  characteristics  of  the 
Leonard- Jones  and  Devonshire  liquid. 


This  work  was  supported  in  part  by  Contract  N7onr -28511  with  the 
Office  of  Naval  Research. 

S 

General  Motors  Corporation  Fellow,  1954-55. 
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The  following  report  simarizes  a method  of  successive  approximations  to 
the  "optimum'’  free  volime  theory  of  liquids.  This  approach  is  reminiscent  of  the 

self-consistent  field  approximation  in  molectilar  quantun  meohinlos.  The  initial 

. 1.2 

approximation  is  the  usual  Lennerd- Jones  and  Devonshire  equation  of  state. 

Detailed  equations  and  procedures  are  set  up  for  determining  the  equation  of  state 

and  the  related  thermodynamic  properties.  It  is  hoped  that  such  calculations  oan 

be  carried  out  in  the  normal  liquid,  as  well  as  the  oritioal  regions. 

Kirkwood  has  suggested  a method  for  Improving  the  Lennard-Jones  and 

Devonshire  "cell  theory"  of  liqviids.  However,  no  oalculations  based  upon  this 

extension  of  the  earlier  theory  have  ever  been  oarried  out.  Qualitatively,  this 

theory  is  not  too  different  from  that  of  Lennard-Jones  and  Devonshire,  in  that 

each  molecule  is  restricted  to  move  within  a cell  bounded  by  other  moleoules  lying 

near  the  centers  of  equivalent  cells  composing  some  simple  crystalline  lattice. 

Although  the  shape  and  orientation  of  the  cells  are  arbitrsury,  we  choose  d''deosdiedral 

oeile  centered  on  the  sites  of  a faco-oentGred-cubio  lattice,  whose  "nearest  neighbor 

Sfeparation"  is  given  by  d^  . 

Let  us  adopt  the  notation  (see  Figure  1); 

(«)  is  tl o vector  from  lattice  site  1 to  lattice  site  j . 

(b)  r,  are  the  displacements  of  the  first  and  jth  molecules,  respectively, 
from  their  sites. 


1.  J.  E.  Lennard-Jones,  A.  F.  Devonshire,  Proo.  Roy.  Soo.  (London)  A163,  53  (1937) 

2.  UTSL.  p.  293 

3.  J.  Q.  Kirkwood,  J.  Cham.  Phys.  18,  380  (I960) 
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•where  all  distances  are  meas\ired  in  imita  of  the  "nearest  neighbor  separation",  d^  . 


For  present  purposes,  , the  pair  interaction  energy  of  moleciiles  1 and 

J , in  units  of  € , is  represented  by  a Lennard-Jones  12-6  potential,  e.g., 

•where  (T  and  € are  adjustable  parameters.  We  define  the  "reduced  temperature" 

T*  , by  T*  - kl/e  . 

In  Kirkwood's  theory  the  "free-Tolume",  , is  given  by, 

(- djr  (2) 

where  *yr  (r^)  is  the  solution  cf  the  integral  equation, 

- c/)  - Eo  (3) 


and 


0(r)-  Z J Vij(ftj)  efrj^dij  (4) 

J>1 

s(r)  - exp  (-■v/r(x’)/T*) 

« exp(-<jf>(r)/T*)  / Jexp(-<^(r)/T*)  dr  (5) 

Eo  *=  J s(r")  <3r 


(6) 
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The  integrations  are  carried  otib  over  cells  of  volume  ^ . In  practice, 

it  is  necessary  to  replace  the  functions,  s(r)  , by  "smoothei"  functions  obtained 
by  averaging  over  all  orientations  of  the  vector  r ; and  the  cell,  A , by  an 
”equi "valent  sphere”.  If,  as  a first  approximation,  we  replace  s(r)  by  a Dirao 
delta  function,  (5(r)  , then  we  obtain  the  results  of  the  Lennard-Jonea  and 
Devonshiro  theory.  Hie  propose  using  this  as  the  first  step  in  an  iterative  prooess, 
to  generate  the  "optimm”  free  volume.  Fortunately,  this  procedure  is  eminently 
suited  to  existing  hign  speed  compublng  machines  such  as  the  QAKAC  at  Wright- 
Patterson  Air  Force  Base. 

When  the  s(r)'s  have  spherical  symmetry  about  the  lattice  sites,  the 
integrations  over  angles  oar  be  performed  exactly.  Thus  we  have. 


- 4,r  z Jstrj)  P ( %rp  rjdri 

i 


J>1  fi 


(7) 


irtiere 


cr 


To  oaloulats  jc^  and  aloe  to  average  <^  (^)  angles,  we  need  the 


result ; 

r .rO 

-j  j 


n 

J*  sin© 

o 


F(Pjj,rj)  d© 


(9) 


f 

-■f^^(rj-^Rij-^r)’-(rj-.Rij-r)*-(rj-  Rjj-r)  } 


for  r<R^j  and  rj<R^j  . 

The  "sums  over  ell  molecules  in  the  system  except  the  first”  can  be  re- 
grouped into  "sums  over  shells  of  neighbors  in  the  lattice".  We  denote  the 
separation  of  the  kth  shell  from  the  lattice  site  of  "the  first  molec\iIe"  by 


I 


15  Ootobar  1954 


, and  the  nuuibei-  of  suok  neighbors  as  . If  -/k  denote  (jr)  as  the  result 

•of  averaging  (_r)  over  all  orientations  rf  , then 


cf>(v)  » ^TT  ^ Yi-jr  fs(r^)  H(a|^;r,ry.)  rj^dr^ 


Ee  “ Arr  JsCr*)  (^)  r^d' 

L 


8tt^  Z s(r)  H(axi  r nTj^)  r rjj^drdr, 


(10) 


(11) 


— 

wnere  ^(^^Ij.)  ®^<1  refer  to  any  mcjacuie  in  the  kth  shell. 

In  the  case  of  the  Ler»nard™Jcnee  and  Devonshire  equation  of  state  it  has 

O O 

proved  useful  to  xntroduoe  the  new  variables  y ' and  y.  “ z*t  . The  trans- 

3 j 

formed  relations  become j 


V2 


s(r)-s(y)-  eyp(-<^<y)/T  y 2n  / exp  (-  cp  iy)/T'*‘)  y dy  , (12) 


H(i 


sy-yj)  = 3yv2[4r'"''’'’y’yj^-  z”r(3;y,yj>]  , 


(13) 


<t"  # > is  the  number  of  raolocules  in  a 


•where  the  " reduced  density  , z 
vol-unie  (^d^)'^  , and 

, r 1/2  Vf  -n  V?  V2,-id  ^ 1/2  1/2  -n  Vg  ^ -rfl 

rCn;y.,y^)  = [(yj^^l-^y  ) ) -(7^- 1-^y  > -t.(yj-i-y)  J(i4) 

The  conbributi.nriF  of  the  flrui;  .-;hell  cf  noighbor.^  for  the  faoe-oentered- 


oubic  lattice  oan  '.hen  be  ■written  ar- . 


i^IsCy,)[^V(a;y,yi) - zV($-,y.yi)J  dy, 
- ;{_s(v)s(yj)[^r('i.y.yj)-zV'r.-/.V,)]<ly^yi 


(15) 

(16) 
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PreviouB  caloiilatloaa  have  shown  that  the  oontributions  to  (y)  and 

by  the  noarest  neighbors  is  far  more  impoirtant  than  for  the  other  shells.  Thus, 
3 Ittle  error  will  result  from  treating  only  the  nearest  neighbors  by  the  more 
detailed  theory  and  using  the  Lennard-Jones  and  Devonshire  for  the  remaining  shells. 
For  a faoe-oentered-oubio  lattice,  the  Lennard-Jones  and  Devonshire  theory 


gives, 

[ ^ -EJlto 

Ce,],„ 


(17) 


(18) 


where 

l(t)  - (1  + 12t  + 26.2t^  + 12t®  + t*)  (1  - t)"^°  - 1 
m(t)  - (1  + t)  (1  - t)"'^  - 1 

are  functions  whloh  have  been  previously  investigated  by  several  workers. 


Thus,  the  oontributions  of  the  2nd  and  3rd  shells  of  neighbors  are  given  by. 


[5Ty)-E.l 


(8.S) 

LTD 


CE.l 


(2i»'5) 


Beoaiise  of  the  short  range  of  the  intermoleoular  forces,  the  effect  of  the 
more  distant  shells  is  oompletely  negligible.  Combining  the  results  obtained 
above  we  have,  ^ 

J S [ff  r (9 ; y ,yj) -zV(z % y,yj)J dy j 

(21) 

+ 96 [7^p(yg*l 


4.  Wentorf,  Buehler,  Hlrsohfelder,  Curtiss,  J.  Chem.  Phys.  1484  (i960) 
5 Kirkwood,  Salsburg,  J.  Chem.  Phys.  21,  2169  (1963) 


TilS-roiE-l?. 

15  Oitober  1964 


r(9;v,yj  )-z^P(5^y,yj)]  dydyj 

■where  s(y)  is  now  the  solution  of  the  modified  integral  equation, 

(i.e,3> 


(22) 


-2-7— ^ (1.2.5)  ^n.a.5> 

(p  (y)  - Eo 


(23) 


Ihe  equation  of  state  of  the  liquid  and  its  internal  energy  per  molecule 
in  units  of  C are  then  given  by. 


pV 


f Xf 

z r A i gyp  V-  <P  (y)/!* ) d ^ ( y)/^2  y Av 

^ ^ y dy 


l-^Ee/SzJ  (24) 


U 


- -m,/^  .T‘/a  .E7ar  . ' 

y*r_A/gyp(-^^)/T*)  d^&)/yr*  v^-^dv  1 
eTTpC-^^/r*')  y’^^oLy  ** 


SCLUTION  OF  THE  INTEGRAL  EQUATION,  (Eq.  23),  BY  AN  TTFRATIVE  PROCEDURE? 

For  a particular  ohoioe  of  the  density,  z , and  the  temperat\n*e,  T*  , 
we  oan  obtain  the  "optimum  fmotion  s(y)"  by  an  iterati-ro  procedure.  For  an 
initial  trial  function,  ^(y)^  s ^be  first  :step  5.n  this  iteration  leads  to  a 
"first  improved"  trial  function,  ^(y)^  . The  next  step,  or  oycl  e,  of  the 
iteration  leads  to  a "second  improved"  trial  function,  s(y)  ; etc.  ... 

■When  the  thermodynamic  properties  calc'jlated  frcm  the  "nth  Impro'ced"  trial 
solution  do  not  differ  soni?J.^;ly  frr»n  -t-hose  calculated  frcm  the  (n-^l)th,  vre 
have  obtained  the  "optimum  function,  s(y)jj  "»  ^ this  sense,  the  procedure 

is  quite  similar  tc  the  "Belf-consi stent,  field  cal c’jJa -cions"  of  atomic  and 


•molecular  quautT-tm  msohanicso 


TTIS-OM-IE 
15  October  1S54 

-7- 


Tho  oyclio  procedure  described  above  oan  be  illustrated  in  the  followir? 

mazuaor  i 
CYCLE  It 

(a)  Choose  a trial  solution  oalled  s(y)^  . 

(b)  fihibstitubing  s(yj)  for  s(y-)  in  Eq.  21  and  performing  the 

Indicated  integration,  vre  oaloulate  the  corresponding  function,  Jo» 


(o)  Substitvrting 


(IAS) 


1 

Jn 


or 


3T— 

Cp(y)  In  Bq.  12 


and  performing  the  indicated  integration,  we  ootain  the  "first  Improved"  function, 

,(X23>i 

i I ff  % 

s(y>] 


‘ ZiT  ferp  (-  I /t'  ) 


CYCLE  2: 


Substituting  for  *(7^)  Eq.  21  and  performing  the 

1,  we  oaloulate  the  corresponding  function  -^1* 


indicated 


(a) 

integration^ 

r 7 y—v  (X£,S)“i  21  ■/' 

(b)  Substituting  [_  (V)  0(y/  Bq.  12 

and  perforaing  the  indicated  integration,  we  obtsdn  the  "2nd  improved"  function. 


iCy^- 


a 


CYCLE  3 1 


rr  /exp(-t3Ty/"-’"^^]a/T’) 

4^ 


CYCLE  n; 


(a) 


(b)  Substituting  L?7y)^^  for 

performing  the  indicated  integration,  we  obtain  the  "nth  inprovod"  function, 

8(yJ  - gyp  <-  [ 1 /t*  ) 


;r-7-ra;i.5) 

<^(y)  in  Bq.  12  and 


ZtT  /exp  (-  [ /t“)  y^^dy 


A 


I 
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This  iterative  procedure  must  be  carried  oxrt  separately  for  each  pair  of  the 
parameters  T and  T*  . Now,  to  obtain  the  thermod;;,caamic  properties  ocrrespondiiig 
to  the  ”kth  improved"  function,  we  must  calculate:  (l)  f u-  ^7  substitu- 

tion  of  ®(y)jj.  f*®*"  s(y)  in  Bq.  22,  and  (ii ) the  derivatives  of  J]^ 

and  arising  in  Eqs.  24  and  25.  If  the  above  iterative  process  is 

oarried  out  sinmltaneoualy  for  several,  choices  of  z and/or  T*  in  the  neighbor- 
hood of  the  values  in  question,  we  can  then  determine  the  derivatives  with  respect 
to  these  peurameters  by  standard  prooediores. 

As  stated  earlier,  wo  shall  choose  s(y)^  " 5 (y)  > "where  (y) 

the  birac  delta  function.  This  choloe  leads  to  the  "first  improved*  function. 


s(y)j  = exp(«l^[z"*  L(y)~  M(y)]  /t*) 

an  GCz.T’*) 


(ec) 


where. 


L(t)  - lit)  -f  -4s  iCt/s) 

128  7k9 

M(t)  ’»  m(t)  ra(t/^)  m U/s) 

.30544 

G(2,T*)-  Jevp(-12.[z^L(y)-2  2^M(v)]/T*)  Y '^dy 

4,5,6 

are  f\mctions  which  have  been  extensively  tabulated  by  several  workers.  • 
Hare,  as  in  subsequent  calculations,  the  integration  interval  is  chosen  tc  be 
(0,  .30544)  o 

The  iteration  from  thi.s  nolnt  on  fellows  the  scheme  outlined  above. 


Unpubl.ished  calculations 
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INTEOBALS ; 

To  o&rry  out  the  proposed  iterative  solution  of  the  integral  equation, 
many  difficult  integrals  must  he  oaloulated  numerioally.  Recently,  we  have 
oalculated  T*)  and  several  related  integrals  such  as  d)G(z,  T*)/^z 

and  2)G(z,  using  the  University  of  Wisconsin  Card  Programmed  Calculator, 

Model  II.  The  integrands,  l(y)  , involved  in  these  oalculabions,  all  had  the 
same  general  characteristics  shown  in  Figure  2. 


Owing  to  this  behavior  cf  the  Integrands,  tTO  adopted  the  following  procedure  for 

the  integrations.  Writing  a typical  in+-evrai  in  the  form 

^944  .0015  >80«  .078  .15C>  .WOG 

J I(y)dy  = fl(v)o.y-  J I(v)dLy  ■*  j I(y>dy  4.  Jl(y>dy+  Jl(y)dy 

o 6 •30S'44  .0015  '078  .1S& 

we,  (a)  asstmed  the  integrand  was  linear  for  the  calculation  of  the  first  two  terms 

(h)  applied  Simpson’s  rule  to  the  last  three  terms  using  the  Simpson’s  interval 
.0015,  .003,  .006,  respeouively. 

It  is  possible  that  a similar  technique  be  applied  to  the  calculation 
of  the  more  oomplicated  integrals  involved  in  the  proposed  iterative  soheone. 
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GRITIQOa  OF  THE  LEKNAIID-JOIJES  ifiiP  DBYONSMRE  THEORY; 

W,  Fiokebt  and  Yvni.  W.  Wood'  ab  L03  Al>imo8  have  rsaently  completed  oalcula- 
tiona  for  Lennard- Jonee  Devonshire  liq’-dd>?  at  very  high  densities  and  temperatures. 

They  fotinl  that  the  Lennard-Jones  and  Devonshire  equation  is  quite  satisfactorj'' 

•under  those  conditions.  Hcrwever,  under  less  severe  conditions,  this  simplified 
theory  is  not  nearly  so  accurate. 

From  experimental  vapor  pressiore  data  on  Xe,  A,  K.r,  N , and  0 we 

have  calculated  the  density  of  a Lennard-Jcnes  and  Devonshire  liquid  in  eqviilibrium 

6 

■sri-th  the  "rapor  phase.  In  performing  these  oalcula-tions  "we  adopted  the  oritioal 
constants  S-j-it  ® .301  and  suggasted  by  Wentorf,  Bue hie r.  Hits chfelder 

and  Ourtiss.^  The  results  of  these  calculations,  together  -with  experimental  ”00- 
exiatenoo  data",  are  presented  in  Figure  Z. 


Ill  IP  -1  ^ 1 ^ T er 

A9*.  2 J*  M J.el. 

Figure  3, 

From  Figure  3 -we  see  that  a Lennard-Jones  and  Devonshire  liqviid  is  more 
.isnso  than  any  of  the  real  liquids  considered.  In  addition,  our  aaloulation* 
indicate  that  the  Lennard-Jones  and  Devonshire  'ralue  for  pV/RT  at  the  critical 


7.  Tables  of  the  Lermard-Jones  ai-.d  Devonshire  Equation  of  state  at  High  Temperatures 
and  Deei3i.*.ioB , YY.  Fi'''.Wa-hc  axid  s7m,  ivocd,  Lcs  iU-smos  Soient’5  fio  Laboratory,  Los 
Aiaiaoa.,New  Mexico  (-anpublishei) 
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poiat  is  .795  compared  to  the  esiperiiaental  ralue  of  about  .292  for  the  rare  gases. 
Ve  also  find  oousiderabls  disagreement  between  experin«  atsQ  and  Leanard-Jones  seid 
DeTonshire  Talues  for  the  -mrlovia  themodynamio  propei'ties. 

CONCLPSIOB; 

In  Tieir  of  the  inadequacies  of  the  Lennard-Jones  and  Deranshire  theory 
mentioned  above,  more  acoxirate  oalculations  in  the  liquid-vapor  oo-ezistenoe 
region  and  near  tlie  oritloal  point  seem  desirable.  Subsequent  oalo\ilations  would 
depend  upon  these  raaialts. 


aft' 


m 


Because  of  our  limited  supply,  you  are  requested  to  return  this  coj^  WHEN  IT  HAS  S.'iJRVE.O 
YOUR  PURPOSE  so  that  it  may  be  made  available  to  other  requesters.  Your  cooperation 

•A*'  * 
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NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINOS.  SPECIFICATIONS  OR  OTHER 

are  Used  for  any  purpose  other  than  in  connection  with  a definttbly  related 

GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  IN<3yjRS 
NO  RESPONSIBILITY,  NOR  ANY  OBUGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  lACENSING  TEE  HOLDER  OR  ANV  OTHER 
RERSON  OR  CORPORATION,  OR  CONVE'mG  ANY  RICS9TS  OR  PERMISSION  TO  MANUFA<J>*UBE 
USE  OR  SELL  AITY  PATENTED  INVENTION  THAT“MAY  IN  ANY  WAY  BE  RELATED  THERETO. 
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